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Abstract 

We calculate the longitudinal (a xx ) and Hall (<r X y) optical conductivities for two-dimensional 
metals with thermally disordered antiferromagnetism using a generalization of an approximation 
introduced by Lee, Rice and Anderson for the self energy. The conductivities are calculated from 
the Kubo formula, with current vertex function treated in a conserving approximation satisfying 
the Ward identity. In order to obtain a finite DC limit, we introduce phenomenologically impurity 
scattering, with relaxation time r. a xx {£l) satisfies the /-sum rule. For the infinitely peaked spin 
correlation function, x(q) oc <5(q — Q), we recover the expressions for the conductivities in the 
mean-field theory of the ordered state. When the spin correlation length £ is large but finite, both 
a xx and a xy show behaviors characteristic of the state with long-range order. The calculation runs 
into difficulty for Q < 1/r. The difficulties are traced to an inaccurate treatment of the very low 
energy density of states within the Lee- Rice- Anderson approximation. The results for c xx (f2) and 
o xy {Q) are qualitatively consistent with data on electron-doped cuprates when SI > 1/r. 



I. INTRODUCTION 



Long-range antiferromagnetic order can have a profound effect on the electronic excitation 
spectrum of metals, opening a gap over some or all of the Fermi surface.- By continuity, 
it seems reasonable to believe that even in the absence of long-range order, finite-range 
correlations may also have an important effect on the electronic excitation spectrum. The 
effects may be expected to be particularly large in two-dimensional systems with Heisenberg 
symmetry, because in this case long ranged order can only exist at temperature T = 0. Even 
in the presence of weak coupling into a third dimension or weak Ising anisotropy, a wide 
range of temperatures will exist where the physics is controlled by the thermally disordered 
magnetic state. Such a state, which following the usual conventions we refer to as a spin- 
density wave (SDW), is believed to occupy a significant portion of the phase diagram of 
electron-doped cuprates.- Extensive experimental studies, including optical conduct ivity,-~- 
Hall effect,— and infrared magnetotransport,- of these materials in the doping range where 
the ground state has long-range SDW order, have revealed signatures characteristic of partial 
gap opening starting at a temperature high compared to the Neel temperature Tpj(x). This 
gap seems to be closely related to that in the SDW state, since the measured quantities 
evolve smoothly across Tn(x).— Furthermore, a recent inelastic neutron scattering study^ 
on the Nd2- a; Ce a ;CuC)4± ( 5 materials found that the spin correlation length £ remains large 
for temperatures high above T/y. This motivates the theoretical proposal that the scattering 
of electrons off thermal spin fluctuations may hold the key to understanding the transport 
data above T N (x).— ^ 

While the experimental phenomenology is clear, our theoretical understanding of this 
regime is incomplete. In a seminal paper, Lee, Rice, and Anderson^ (LRA) proposed a 
model for the study of electron dynamics in the presence of long but finite ranged density 
wave order. In this model, electrons are coupled to quasi-static (relevant frequencies less than 
ksT) order-parameter fluctuations, resulting in suppression of the single-particle density of 
states at low energies, a phenomenon sometimes referred to as a "pseudogap" . The Lee- Rice- 
Anderson analysis was generalized by Sadovskii^ and then was extended to two-dimensional 
systems close to the antiferromagnetic instability by Vilk, Tremblay, and co-workers.— ^ who 
argued that such long but finite ranged antiferromagnetic fluctuations controlled important 
aspects of the physics of the electron-doped cuprates. In a further theoretical development, 
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Schmalian et al argued that the electron Green's function can be exactly calculated for 
the two-dimensional LRA model of electrons with a cuprate band dispersion scattered from 
antiferromagnetic spin fluctuations,— generalizing the method first used by Sadovskii in 
the study of one- dimensional charge-density fluctuations.— Tchernyshyov^ analyzed the 
underlying assumptions of the Sadovskii's solution, and argued that it should be used with 
caution in the generical two-dimensional situation. However, he concluded that in particular 
regions of momentum space termed "hot spots" , the method could be safely applied, and it 
is near these momentum points that Schmalian et al found pseudogap behavior in agreement 
with the previous work of Vilk and Tremblay.— ^ There has been an attempt to calculate 
the conductivity using this method.— However, the restriction to the vicinity of the hot 
spots makes the Sadovskii solution unsuitable for the study of transport properties in two- 
dimensional systems, because a summation over the entire Brillouin zone is needed. A 
generalization of the LRA theory to transport phenomena is required. 

In this paper, we provide the missing generalization. We use the two-dimensional LRA 
model in which electrons are coupled to themal (quasi-static) antiferromagnetic spin fluctua- 



tions to study the optica 



above T N . As in Refs 



and Hall conductivities of electron-doped cuprates at temperatures 



10Ml2| , we calculate the electron self-energy in the leading order of per- 



turbation theory. The new feature of our work is a calculation of the current vertex function 
in a conserving approximation.—^ We find that although the vertex function corresponding 
to the LRA self energy leads to a conductivity which fulfills the /-sum rule, the dynamic 
(H — ) current vertex function has unphysical features at low frequencies; leading in some 
cases to an unphysical negative conductivity in the very low frequency region. The difficulty 
is traced to an incorrect treatment of the subgap density of states in the LRA calculation. 
We discuss ways of curing the difficulty and also present results at higher frequencies which 
are not significantly affected by the problem. 

The rest of the paper is organized as follows. In Sec. [Til we use the spin-fermion model to 
motivate the LRA model, calculate the electron self-energy in the leading-order perturbation 
theory, and discuss the pseudogap phenomenon in the resulting single-particle spectral func- 
tion. In Sec. PU we study the optical conductivity with a proper treatment of the current 
vertex function. In Sec. IIV| we study the Hall conductivity, developing a calculation scheme 
which can reproduce the mean-field result in the proper limit. In Sec. [V] we summarize 
our results, and discuss the implications. Some technical details and a brief summary of the 
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mean-field theory can be found in various Appendices. 
II. MODEL AND ELECTRON SELF-ENERGY 
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FIG. 1: Quadrant of two-dimensional Brillouin zone showing the Fermi surface (solid line) for 
electron-doped cuprates, its translation ('backfolding') by wavevector Q = (7r,7r) (dashed line), 
'hot spot' p* and other Fermi surface points referred to in subsequent discussions. 

In this section, we present the model, the basic approximation we use, and results for 
the self energy. The results reproduce those previously derived^ - — and are presented 
here to establish notation. Our starting point is electrons moving with a two-dimensional 
dispersion chosen, for concreteness, to resemble that believed to be relevant to high-T c 
superconductors:— 



where t = 0.38eV, t' = 0.32t, t" = 0.5t', and /x is the chemical potential. Figure [T] shows the 
resulting Fermi surface for electron doping x ~ 0.16 (solid line). Throughout this paper, we 
choose units such that the lattice constant a — 1, H — 1, and measure energy in units of t, 
unless otherwise stated. 

We represent the effects of magnetism via the spin-fermion model, which has been used 
extensively in the study of itinerant electrons in systems close to or in long-range magnet- 



e p = —2t(cosp x + cospy) + At' cosp x cosp y — 2t"(cos 2p x + cos 2p y ) — fi, 



(1) 
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ically ordered states.— — It is a low-energy effective theory with an intrinsic cutoff energy 
A, and is conveniently formulated as an effective action^ 

rP rP 



where c pa is the fermionic field operator, Gq (p, r) is the inverse of the bare fermionic Green's 
function, S q is an emergent field describing collective antiferromagnetic spin fluctuations 
which are important to the low-energy physics, Xo(q>^) — Xo/[£~ 2 + (q — Q) 2 — ( w A's) 2 ] 
is the bare susceptibility in the spin-fermion model with Q = (tt,tv), s q = ^2 p c' p+qa (T a pCpp 
is the spin density operator of electrons with cr the Pauli matrices, and g is the effective 
coupling constant between electrons and spin fluctuations. 
Go(p, t) is the Fourier transform of 



where e n = (2n + 1)ttT, and we have explicitly introduced the impurity scattering rate l/2r, 
which will be discussed in the next section. Eq ([3]) has been extensively used in studies 
of the fluctuation conductivity close to the Peierls transitional^ and the superconducting 
transition.—^ When two different scattering processes, spin fluctuations and impurities, 
are present, it is necessary to consider their interference.— 1 ^ The renormalization of the 
spin-fermion interaction vertex g by impurity scattering and that of the impurity scattering 
relaxation time r by the spin-fermion interaction are discussed in Appendix [A], where it is 
demonstrated that both renormalizations are finite. Thus, as long as we keep g and r as 
adjustable parameters of the calculation, we can neglect their mutual renormalizations. 

The third term in Eq (j2J) represents the interaction between electrons and spin fluctua- 
tions, and effectively arises from an interaction Hamiltonian, 



p.q 

In this paper, we consider the state without long-range order, < S q >= 0, and assume that 
the spin fluctuations are isotropic. As a result, the spin indices on electrons are irrelevant 






(4) 
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for the calculation of charge transport coefficients. After properly redefining g to account 
for the three S directions and two electron spin projections, the interaction Hamiltonian can 
be written as 

which bears the form of the electron-phonon interaction with S q playing the role of the 
phonon field operator A q = a q + aL q .— There are important differences between the spin- 
fermion model in Eq (j3J) and the electron-phonon problem. The phonon degrees of freedom 
are extrinsic to electrons. Because of the large mass of nuclei compared to electrons, the 
Migdal's theorem applies, which allows the electron-phonon interaction in conventional met- 
als to be treated in a controlled manner.— However, the spin-fluctuation degrees of freedom 
in the spin-fermion model are intrinsic to electrons. For the Migdal's theorem to be applica- 
ble, one usually resorts to one or another variant of the large- N limit where N is the number 
of fermion flavors^ or the number of hot spots.— As in the electron-phonon problem, the 
bare spin- fluctuation propagator xoio.^) is renormalized by creation and annihilation of 
electron- hole pairs, which leads to the Landau damping term, iu/u s {. In the random phase 
approximation, the renormalized spin-fluctuation propagator has the form 

s + (q - Q) + K^/Wrf 

where cu s f/£ 2 sets the energy scale for spin fluctuations, and can be expressed as combina- 



tions of the parameters in Eq (J2]) (see e.g. Ref |24j). We note that Eq (jHJ) has the same 
form as that proposed phenomenologically by Millis et a/,— can be obtained from the self- 
consistent renormalization theory,—^ and has the generic form in the theory of quantum 
phase transitions involving itinerant electrons.— i22j2& We thus argue that the applicability 
of Eq (jHJ) is independent of microscopic details and approximations involved in deriving it. 
The remaining question is to calculate effects of the interaction in Eq (jSJ) on fermions. 

When the temperature T is relatively large compared to u s f/C, 2 , the dynamic term in Eq 
(JBJ) can be neglected (for more discussion, see Appendix |Bl and for a related discussion in the 



context of superconducting fluctuations, see Ref |29|). In terms of Matsubara frequencies, 
the static spin-fluctuation propagator is written as 

X(q, iu n ) = f° Ty^^O, (7) 

C + (q- Q) 

which is the two-dimensional generalization of the LRA model. 
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FIG. 2: Feynman diagrams used in this paper, (a) The Lee-Rice-Anderson approximation to the 
electron self energy, (b) the current vertex function corresponding to the Lee-Rice-Anderson self 
energy and (c) diagrams needed for the longitudinal conductivity. The thin solid lines in (a) and 
(b) represent the bare electron Green's function Go in Eq ([3]), the thick solid lines in (b, c) represent 
the dressed Green's function G in Eq (|16|) . the wavy lines are the spin- fluctuation propagator x 
in Eq ([7]), and the solid circle in (c) represents the current vertex function T J and is calculated 
according to (b). 

Using Eq ([3]) and Eq ([7]), we calculate the electron self energy to leading order in g, shown 
in Figure [2] (a), 

2 ^ [ dq xo 1 



£(p,ze„) — g T 



(2tt) 2 (q - Q) 2 + £" 2 ie n - e p+q + *sgne n /(2r) 



(8) 



As discussed above, there is a cut-off energy scale A below which the spin-fermion model is 
defined. Since the above integral is convergent at large q, we assume that all energies under 
consideration are below A. To proceed, we change q — > q+Q and write q = (q\\, q±) where q\\ 
and q± are the components parallel and perpendicular to v p+ q, respectively. The gy-integral 
can be done by the residue method, and the remaining g^-integral is elementary. The 
retarded self-energy is obtained by analytical continuation to the real axis via ie n — > u + i5 
with 5 a positive infinitesimal, 



P,u J ) = —f( UJ ~ £p+Q ' 



(9) 



Here we introduce the effective coupling constant A and energy scale Eq defined by 

A = g 2 T X o/(2ir), (10) 

£0 = Vp+gC 1 , (11) 



where the weak p-dependence of eo will be neglected in subsequent calculations. The function 
F is given by 



F{z) = / dx * = - 1 In 1 + ^ + V gH . (12) 







In the limit \ jr — > 0, we reproduce the result of Refs lOlJllj ] 



E^-(p, u) = =|sgn(cu — s p +q) In ° = 

V(w - e p+Q ) 2 + eg L v/(o;-£ p+ q) 2 + ^- |w-£ p+ q| ^ 



-17T/2|. 

We note that Im£^(p, a;) < as expected and have verified that ReS^j.(p, u) and 
ImE^(p, oS) are related by the Kramers- Kronig relation, 

frfljKp,-,) = lp f^' mS " (P ^>. (14) 

Using the one-dimensional analogue of Eq ([8]), one obtains the result of Lee et al,— 

= f. /9 ^. , (15) 

w - e p+Q + ?/2r + 2£ 

where A 2 = g 2 Txo£,/2. We note that Sf D (p, w) has a simple pole at u = s p+ q — i/2r — ieo 
in the lower-half w-plane, as expected. 

We now discuss the physical content of the results. The retarded Green's function is 

G *(p> ^ = 7^ VR( v (16) 

us — e p + i/2r — L-"(p, oj) 

and the spectral function A is 

A(p, u) = -2ImG K (p, u). (17) 

For comparison to previous results, we will sometimes present results based on Eq ( 1T5|) 
(i.e. for the model without an explicit additional impurity scattering). We denote the 
corresponding Green's function and spectral function as Gf* and Aif, respectively. 

The situation is particularly simple in the one-dimensional case. If we linearize the 
dispersion near the Fermi level e p — > vp, e p +q — » —vp, measure momenta relative to the 
Fermi momentum ±kp, and assume Q = 2k f, then 
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exhibiting a gap of size A 2 + (vp) 2 broadened by the impurity scattering rate and by the 
finite correlation length (parametrized by Eq). In obtaining this result, it is crucial to use 
the bare Green's function in Eq OH]). Self-consistent one- loop approximations (and related 
approximations such as the fluctuation-exchange approximation (FLEX)) do not obtain a 
pseudogap. 

In the two-dimensional case of main interest here, the situation is more complicated 
because the Green's function depends both on position on the Fermi surface and on dis- 
placement of the momentum away from it. However, a few general statements can be made. 
We note that Eq ffT5]) can be written as E^(p,cj) = £ [(\/i£l)F(uj/i£ — e p+ q/ieq)\. Thus, 
the spectral function shows scaling behavior: £ Aif(p, UJ / e o) is invariant if energies and 
frequencies are measured in units of £ at fixed p and X/sq- 

Precise results can be obtained in the limit that the spin- fluctuation propagator x(q) is 
infinitely peaked at Q (the Kampf-Schrieffer model^) 

XKs(q, iu n ) = xo<$n,o<Kq - Q). (19) 

Eq OH]) gives 

££ 5 (p, u) = A ' , (20) 

u - £ P+Q + i/2t 

where A 2 S = g 2 Txo/ (27r) 2 . Substituting this self-energy into Eq ([TBI) , one obtains the Green's 
function in the mean-field theory of the SDW state (the diagonal matrix elements in Eq 
( 1C3D ) without introducing a condensate.— In the next two sections, we shall extend this 
conclusion to the optical and Hall conductivities; the mean-field expressions for <j xx and o xy 
can be obtained from Eq (fT9l) in the leading-order perturbation theory. 

An important role in subsequent discussions is played by the "hot spots", momenta p* 
such that both p* and p* + Q are on the Fermi surface (Figured]). At these points, the 
density of states is most strongly reduced from the non-interacting value. The structure of 
the spectral function at the hot spots is parameterized by a gap scale A pg and the scaling 
arguments of the previous paragraph show that A pg /£o depends only on A/^q. In the limit 
A/£q 3> 1, A pg /£ ^> 1, and is determined by the equation 

\ g = Re£^(p*, A p g). (21) 
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FIG. 3: (a): comparison of the spectral functions A(p* ,uj) (solid line) and Aif(p*,uj) (dashed line) 
at one hot-spot, p* (see Fig. [TJ. The parameters are £o = 0.2i, A = 2.5£q, and 1/2t = 0.05i. (b): 
the spectral function Aif(p,u;) at p* (solid line, black online), p c = (tt/2,it/2) (dashed line, green 
online), pi = (0.43-7T, vr/2) (dash-dotted line, blue online), and p2 = (0.57-tt, it/2) (dotted line, red 
online). The parameters are £o = O.lt and A = 10£q. 

To leading logarithmic accuracy, we find 



with a (weak) logarithmic dependence on e . This equation should be contrasted to the 
one-dimensional result A 2 = g 2 Txo£,/2- 

Panel (a) of Figure [3] shows the spectral functions A(p*,u) (including impurity scattering, 
solid line) and Aif(p*,u) (no impurity scattering, dashed line) at the hot spot p*, for 
l/2r = 0.05t, Eq = 0.2t, and A = 2.5£q. Both curves show suppression of the spectral weight 
at low frequencies. We define the pseudogap A pg as half the distance between the two peaks 
on the corresponding curve and see that the two curves have roughly equal pseudogap values, 
~ OAt, slightly larger than that predicted by the asymptotic result in Eq (122|) . m 0.351 
This panel thus demonstrates that we can use either Eq (Q or Eq ([TBI in discussions of the 
pseudogap in the single-particle spectral function if impurity scattering is reasonably weak. 

Panel (b) of Figure [3] shows Aif(p,ui) for several p (see Fig. [Q, including a hot spot 
(solid trace, black online), a point p c which would be at the center of the "hole pocket" in 
the SDW state (dashed trace, green online), a momentum pi far from the hot spot but on 
the noninteracting Fermi surface and a momentum p 2 which would be near the back side of 




1/2 



(22) 



Thus in the two-dimensional case, in the limit X/eq ^> 1, A pg is determined mainly by y/\ 
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e =0.1t 
e =0.2t 
e =0.3t 



FIG. 4: The spectral function ^4(p*,o;) at the hot spot p* for 1/r = O.lt and fi = 0.175t. (a): 
X/t 2 = 0.1 and different values of (b): X/t 2 = 0.3 and different values of Eq. 

the hole pocket. At the hot spot, one observes two peaks, symmetrically disposed around 
the chemical potential. At the center of the "hole-pocket", one also sees two identical peaks, 
but this time not centered at the chemical potential. At the other two momenta, one sees a 
large peak indicative of a conventional Fermi liquid quasiparticle and a small 'shadow peak' 
at the location of the other quasiparticle state. All of these features may be understood in 
terms of a broadening of the mean-field solution. 

In Figure H] we investigate the dependence of the spectral function on the parameters A 
and Eq. Each panel shows the spectral function at the hot spot, computed for a fixed A and 
several different Eq. The gap scale (defined from the peak to peak distance at the smallest 
Eq) increases with increasing A. As Eq is increased, the low-energy density of states increases 
(gap fills in) and at larger Eq, the gap magnitude (defined from the peak separation, when 
visible) decreases, but at a rate slower than the increase of the low-energy density of states. 
Thus the suppression of the pseudogap has more to do with the gap filling in than with a 
gap decrease. 

Comparing our results to data suggests that it is reasonable to associate an increase in 
temperature with an increase in the parameter Eq (i.e. a decrease in the correlation length 
£ as observed in Ref j^]), while the increase of A produces effects similar to those observed 
in electron-doped cuprates when doping is decreased. To qualitatively relate theory to 
experiment, we therefore fix the chemical potential and model changes in doping by changes 
in A and changes in temperature by changes in Eq. 
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III. CURRENT VERTEX FUNCTION AND FREQUENCY-DEPENDENT LON- 
GITUDINAL CONDUCTIVITY 



The longitudinal conductivity a xx is given in terms of the polarization function II as 

u P (in n ) + u D 

a xx {iU n ) = (23) 

where the paramagnetic P and diamagnetic D contributions to the polarization function are 
given in terms of the current vertex T J as 

^-pWpG(p, ie n )T J x (p, p + q, ie n , ie n + itt n )G{p + q, ie n + ift„), 

(24) 

where gq = e 2 /h is the conductance quantum and v* = de p /dp x , and 

U D = IoqT £ J j^y 2 e; x G(p, iO, (25) 

where et x = d 2 e p /dp x . 

The magnitude of the current vertex function T J is related to the relative sizes of the 
frequency and momentum dependence of the self energy. We have seen in the previous 
section that the momentum dependence of the self energy is not negligible and thus expect 
the current vertex correction to be important. An important constraint on calculations is 
the Ward identity following from current conservation; this ensures that the conductivity 
obeys the "/-sum" rule. The Ward identity relates the density vertex T p and current vertices 
T J to the electron propagator via 

G _1 (p+q, ie+iQ) -G -1 (p, ie) = i£lT p (p, p+q, ie, ie+iQ) - q-f J (p, p+q, ie, ie+iQ). (26) 

Taking the q — > limit with Q fixed to gives 

hm f J (p, p + q, ie, ie) = v p + dE & l6 \ ( 27 ) 
9^0 F dp 



To obtain the vertex function, we follow the procedure outlined in Refs |41|J42|: insert the 
free vertex on each bare electron line in the diagrammatic expansion of the electron Green's 
function, and then amputate the resulting diagrams. The diagrammatic expansion for the 
current vertex function is shown in Fig|2] (b), and the corresponding analytic expression is 
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(henceforth we drop the superscript J and remove one of the two momentum arguments 
because we deal only with the current vertex in the q — > limit) 



r x (p, ie n , ie n + iQ n ) = v x p + g 2 T J 



Xo 

(27r)H- 2 + (q-Q) 2 

x C (p + q, ie n )G (p + q, ie n + iVL n )v p+q . (28) 



Approximating v p+q ~ v p+ q, Eq ( 128|) is evaluated as 

r / • • , -n x x, S(p,ze n ) -S(p,ze n + zfi n ) x 

r.(p, > e „ + rfi.) = « P + lfi „ + i[sgn(e „ + fi „ ) _ sgM , J/2 /V + o. (29) 

which is consistent with Eq (127)1 because Eq (JSJ) shows that dTi(p,ie n )/dp a = 
— fp + Q<9E(p, itn)/d(ie n ). We therefore conclude that Eq (!2"9j) for the current vertex function 
is a conserving approximation. 

In the study of the fluctuation conductivity near the superconducting transition, the 



Aslamazov-Larkin (AL) contribution is usually important (see e.g. Ref 29|)- Its calculation 
requires the inclusion of the dynamic term in the spin-fluctuation propagator x(q, u). We 
calculated the AL contribution using Eq (JSJ), and found that it is negligible, because x(q, k>) 
is peaked at a finite momentum Q. 

The physics of the vertex correction may be understood by comparison to the mean-field 
solution in the ordered state. To demonstrate the main issues with a minimum of notational 
complexity, we discuss the one-dimensional model, in which we linearize the dispersion about 
the Fermi energy, measure momenta from the Fermi momentum, and assume the ordering 
wave vector Q = 2kp. The mean-field solution is characterized by normal (Gmf ~< c p c p >) 
and anomalous (F ~< c p c p+ q >) Green's functions given for right (a = +) and left (a = — ) 
moving electrons by 

G a MF (p,iu) = — * 30 

ur + v z p z + A z 

F a MF {p,iuj) = -— —A—-, (31) 
and the conductivity is given schematically by 

a oc ^Tr [GG - FF] , (32) 
with the trace over frequency, momentum and left/right index a. 
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Turning now to the theory in the fluctuation regime, we have (for u> > 0) 

G L ra{p,iu) = 5 5-5 77: 771 n F7 \T' ( 33 ) 

L0 Z + u z p z + — l[[lU) — vp){7^ + Eq) + 27(«w + vp)\ 

The absence of long-range order means that F = 0, so that cr is evaluated directly from Eq 
(J23l) while use of Eq ( |29l) for the vertex function and Eq ( fl5l) gives 



- U ( 1 - 7 71 ^777 71 77 I , (34) 

(iw + ifi + vp + i + £0)) (iu + vp + i(± + e )) 1 



if sgn(a; + Q) — sgn(a;) = +, and 



A 2 (l + 



2e 



(iw + iti + vp + z + e )) (iw + vp - i + e )) 

if sgn(ci; + = + but sgn(a;) = — . 

Substituting into Eq ([25]) . we see that the first of the two terms in the vertex function 
reproduces the GG term. The second of the two terms reproduces the FF contribution, 
which, in the ordered state, carries the coherence factors which for example distinguish 
antiferromagnetism from superconductivity. Thus the vertex correction does what is required 
to produce the correct form of the near gap conductivity. However, we see that in addition, in 
the physically crucial sgn(o;+fi) 7^ sgn(w) regime there is an extra term, of order e / t) 
which diverges asfi + l/V— y but is unimportant for \Q + l/r| > Eq. The structure of this 
term is a defect of the Lee-Rice- Anderson approximation. We believe it occurs because this 
theory produces an incorrect form for the subgap density of states, which should vanish as 
frequency 00 — > 0. Indeed in the one-dimensional case it is known that the low-frequency 
density of states is due to amplitude singularities in the flucuating order parameter, which 
become exponentially rare at low frequencies.— 

The problem can also be cured by a self-consistent treatment such as FLEX, but this is 
known to give an incorrect form for the pseudogap density of states.— We have not been able 
to identify a consistent and physically reasonable cure for the divergence which is applicable 
also in two dimensions, so we adopt the expedient of introducing an impurity scattering 
which cuts off the divergence. We shall see, however, that the theory can still produce an 
unphysical dip in the low-frequency conductivity. 

Returning to the two-dimensional model of primary interest in this paper, we combine 
Eqs. ( l23ll24ll25lT29l) . perform the analytical continuation, and obtain 

Rea xx (n) = Rea£(n) + Rea^fi) + Rea£")(fi), (36) 
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where 



Reagan) =a Q 



Q 2 + 1/r 2 J (2tt) 2 p p+q 

duf(u) - f(u + n) 



x 



[A(p,w)+A(p,w + n)], (38) 



and 



Reai{ n \n) = 2a Q 



1/r 2 f dp 



W + 1/t 2 J (2^p p P+Q 

duj f(u) - f(u + Q) ReG(p, u + Sl)- ReG(p, g) 



(39) 



where f(x) is the Fermi function. In the calculation, we assume that the most important 
effect of temperature is on £ (or Eq), and neglect thermal broadening of the Fermi function. 
As a result, f(x < 0) = 1 and f{x > 0) = 0. In the limit r — » oo, Reuii remains finite, 
Rea£ 7) = AStf(ft), and Rea£ /7) (ft) -»■ 0. Thus, Re(x£ 7) is the divergence discussed above. 
Its weight AS* is given by 

We find that AS < for the band dispersion appropriate to cuprates, since v*Vp + Q < in 
most part of the Brillouin zone where A(p, 0) is appreciable. 

Figure [5] shows Kea xx (Q) calculated from Eq fl36|) with 1/r = O.lt and \x = 0.175t for 
different values of A and e - All curves in Figure [5] show anomalous low-frequency behavior, 
arising from Eq ( 1381) . Since AS < 0, the DC limit can be made negative (not shown 
here) for larger values of A or r. As Q increases, this anomalous contribution is quickly 
suppressed due to the prefactor Q2+i/ T a ■ I n Figure [5], we see that Recx x . x .(f2) behaves as 
expected for Q > 3/r, and we shall concentrate on this regime. In this regime, for small Eq 
(eo/t = 0.05,0.1,0.2 in panel (a) and Eo/t = 0.1,0.2,0.3,0.4 in panel (b)), there are peaks 
around 2A pg as determined from Figure HI This peak structure is reminiscent of that in 
mean- field calculations, as shown in Figure ITOl (a). The peak becomes weaker for larger Eo 
(smaller £). We now use the association of A and e with x and T as discussed in Sec. [TT] 
to relate these results to experimental observations. Comparing the two panels in Figure [5] 
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FIG. 5: The longitudinal conductivity per plane evaluated from Eq (I36j) with 1/r = O.lt and 
fi = 0.175i. (a): X/t 2 = 0.1 and different eo- (b): X/t 2 = 0.3 and different eo- To convert to 
physical units, one must multiply the calculated result by the conductance quantum, oq = e 2 /h, 
and divide it by the inter-plane distance. 

suggests that at low temperatures, Rea xx has an optical peak, the peak position decreases 
with doping, the peak vanishes at some temperature T*, T* increases with underdoping, 
and for fixed doping, there is a spectral weight transfer from high-frequency region to low- 
frequency region, as T is increased. 

We have verified numerically that the calculated conductivity obeys the /-sum rule 

/ dnRea xx (Q) = -U D . (41) 
Jo ^ 

To see this analytically, we consider the case of small 1/r, such that KeaxJ 1 ^ can be 
neglected and ^RecT^^O) can be approximated as 

r -Rea^m « AS/tt. (42) 

At the same time, J*^° — Recrii can be evaluated using the Kramers-Kronig relation 
between A(p,u) and ReG(p,cj), and the result is 

/°° (70 I" rlrt f°° r\f 

^ — Re°i2(n) = J J^y v >P - v p + q) J ^ -m^P, ^G(p, e). (43) 

ITd is found to be equal to the sum of the above two equations. Since R,ea xx (Q) is an even 
function of Q, we obtain Eq (jUJ). 
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The conductivity a xx (fl) in the mean-field theory of the SDW state (Eq fjC4j) ) can be ob- 
tained in the present leading-order perturbation theory, by simply substituting Eq (120]) into 
various Green's functions in Eq ( I36p . The vertex corrections are crucial in this derivation; if 
we had neglected the vertex corrections, we would effectively have neglected the off-diagonal 
terms in Eq (]C3p (the terms proportional to A). The conductivity in the mean- field SDW 
state for A = 0.3t (solid line) and A = 0.6t (dashed line) is shown in Fig [10] (a). The low- 
frequency parts of these curves are well described by a Drude peak without any anomalous 
dip. One possible reason is that in the limit r — > oo, both A and G R in the mean-field theory 
are singular, unlike Eqs (0, [TB]) which are finite due to scattering from spin fluctuations. 



IV. HALL CONDUCTIVITY IN THE PERTURBATION THEORY 



In this section, we develop the formalism for calculating the Hall conductivity a xy in 
the leading order perturbation theory. Th e calculation of a xy in the self-consistent Born 



approbation can be found in Refs WW At the level of approximation emp.oyed here, 
we find it easier to apply the method developed in Ref [47| to the conductivity diagrams 
shown in Fig|2](c). Rewriting the diagrams in Fig|5](c) in terms of the bare Green's functions, 
replacing every electron momentum p in the loop, according to the minimal coupling rule, 
by p— ^A, and expanding the resulting diagrams to first order in A, we obtain the diagrams 
shown in Figure |6] in which the intersections where the magnetic field lines denoted by B 
meet the dressed Green's functions G(p, ie n ) represented by thick solid lines are the dressed 
magnetic vertices, which are calculated according to Fig|2](b). 

Summing all the diagrams in Figure M and expanding the resulting expression up to first 
order in k, we find that (1) the terms independent of k vanish, and (2) the terms of first 
order in k depends on B = ik x A = Bz, signaling gauge invariance. The Hall conductivity 
on the Matsubara axis can be expressed as 

<j xy (iQ n ) = \°Q^^f n J ^{Sx -S 2 + S 3 + S 4 - S 5 }, (44) 
where $o = ithc/e is the superconducting flux quantum, 

5 1 = TY J {d y Vpi - d y FX) {GdM+) - Gi+W), (45) 

ic,, 

5 2 = TY J {d^y p - d * T >p) {Gd y S(+) - Q(+)d y g), (46) 
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FIG. 6: Diagrams for the Hall conductivity. Thick solid lines represent the dressed Green's function 
G(p,ie n ), thin solid lines represent the bare Green's function Go{p,ie n ), and wavy lines represent 
the spin- fluctuation propagator. The vertices with one dashed line are associated with = 
dEp/dp a , the vertices with two dashed lines are = d 2 e p /dp a dp/3, and the vertices with three 
dashed lines are e^ 7 = d 3 e p /dp a dp/3dp 1 . The Greek indices, a, ■ ■ ■ , refer to directions of the 
external fields, E, J, and B. 



Ss = TJ2( T > V P ~ r X) (9 x S(+)d y g - «+)), (47) 
Si = 9 2 T 2 J2v y P GG(+) J 7^x(p-q)(<X-^X)(^o(+)^-^o(+) 2 ), (48) 

S 5 = g 2 T*Y,v x p gg(+) j ^x(p-q)(«-«)^(+)5o 2 -^o(+) 2 ). (49) 



and 



In writing these equations, we have used short-hand notations: Q = G(p,ie n ), Q(+) = 
G(p, ie n + itt n ), Tp = F x (p, ie n , ie n + iQ n ), £ = G (q, ie n ), and G (+) = G Q (q, ie n + iVt n ). 
Using the spin-fluctuation propagator in the Kampf-Schrieffer model, Eq fll9p . after a lengthy 
calculation, we can show that Eq (1441) reduces to that in the mean-field theory of the SDW 
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FIG. 7: The Hall conductivity per plane divided by B, Ima xy {Vl) / B , calculated from Eq (|D1|) 
with 1/t = O.lt and n = 0.175t. (a): X/t 2 = 0.1 and different values of e . (b): X/t 2 = 0.3 and 
different values of £o- To convert to physical units, one must multiply the calculated result by the 
conductance quantum and the in-plane unit cell area, and divide it by the superconducting flux 
quantum, <3>o = hc/2e, and the inter-plane distance. 

state (Appendix IC]) .— The proper treatment of the vertex functions as discussed here is 
crucial in arriving at this conclusion. 

The frequency summation in Eq (144j) is standard.— The physical observable a xy (Q) is 
obtained by analytical continuation iQ n — > Q + i5. The expression for lma xy (Q) is quite 
cumbersome. Here, we focus on the results, leaving detailed expressions to Appendix iDl 

Figures [7] and |8] show Ima xy (Q) calculated from Eq f ]Dlj) for 1/t = O.lt, /i = 0.175t, 
and various values of A and Eq. The frequency scales in Figures [7] and |8] are selected to 
highlight the frequency window most relevant to experiments^ where a xy was measured for 
Q > 0.3t 1000cm -1 . The low-frequency part ( Q < 1/r) of lma xy (fl) suffers from the 
same difficulty as does Rea xx (Q) due to the perturbative nature of the calculation and we 
do not show results in this region. 



Comparing to Ref 
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suggests that the present calculation captures important features 
of data. First, Imcr^ can be made negative at low frequencies, although we start with 
a single-band model with a hole-like Fermi surface. A previous study suggested that the 
appearance of the negative lma xy in electron-doped cuprates is a signature of the long-range 
spin-density wave order.— Our findings here suggest that fluctuating order can also explain 
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(a) (b) 




£2/t QJt 

FIG. 8: The Hall conductivity per plane divided by B, Ima xy (Q)/B, calculated from Eq (|Dip with 
1/r = 0. li and (J, = 0.175t. (a): £o/t = 0.1. (b): £o/t = 0.3. To convert to physical units, see 
Figure [7J 

this behavior. The vertex corrections shown in Figure [6] are important for this conclusion. 
We found that lma xy remains positive in the entire frequency range, if only the diagrams 1-4, 
11, and 12 in Figure [6] are kept and the vertex corrections to magnetic vertices are neglected. 
Kontani and co-workers, using the FLEX approximation, also emphasized the importance 
of the magnetic field vertex corrections.— As shown in Figure [7J lma xy (u) has relatively 
sharp peaks for small e (e /t = 0.05,0.1 in panel (a), and e /t = 0.1,0.2 in panel (b)) 
around 2A pg , showing precursor effect to that obtained from a mean- field calculation shown 
in Figure fTUI (b). The peak gradually vanishes as Eo increases. Furthermore, fixing A (or 
doping x) and increasing Eq (or temperature T), lma xy (Q) increases from negative to positive 
at low frequencies and decreases at high frequencies. This is qualitatively consistent with the 
trend observed experimentally in electron-doped cuprates in the underdoped regime.— From 
Figure El we see that at fixed Eq, Ima xy (Q) decreases with increasing A (or decreasing x) from 
positive to negative at low frequencies, and increases with increasing A at high frequencies, 
again qualitatively consistent with data.— However, in our study, lma xy (Q) remains positive 
at high frequencies, inconsistent with data.— More quantitatively, if (following the discussion 
of the longitudinal conductivity above) we assume that X/t 2 = E /t = 0.1 is a reasonable 
representation of cuprates at 0.12 electron doping, we see that the predicted zero crossing 
in a xy occurs at Q ~ 0.15 — 0.2eV, again semiquantitatively consistent with data. However, 
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our calculation exhibits more temperature dependence than is found in data. 
V. SUMMARY 

In this paper, we used the Lee-Rice-Anderson model to study two-dimensional electrons 
scattered from static antiferromagnetic spin fluctuations, with potential applications to 
electron-doped cuprates in the underdoped regime where the long-range spin-density wave 
ground state is expected. Our theory is in a sense complementary to that of Kontani et al 
who used a fluctuation-exchange approximation most applicable in the overdoped region.—^ 
The theory has two important parameters: A which controls the gap amplitude, and Eq which 
represents the effect of non-zero temperature. 

We first discussed single-particle properties with the self-energy calculated in the leading- 
order perturbation theory. There is pseudogap opening for relatively small Eq which is related 
to the spin correlation length £ via Eq = vp/C, where vp is the Fermi velocity. As Eq increases, 
the pseudogap is gradually filled in with a moderate change in the size of the pseudogap. The 
value of the pseudogap is primarily determined by the coupling constant A between electrons 
and spin fluctuations. We assume that A is increasing as underdoping, and Eq is increasing 
as increasing the temperature. The conductivity Kea xx (Q) is calculated in a conserving 
approximation which respects the /-sum rule. The current vertex has unphysical low-energy 
features. We found that in order to obtain a finite DC conductivity, it is necessary to include 
impurity scattering. However, even with impurities, the low-frequency part of Rea xx (Q) still 
behaves anomalously. As frequency Q increases larger than 1/r, the anomalous contribution 
is quickly suppressed. Rea xx (Q) is characterized by a peak around twice the pseudogap value 
for relatively small Eq (large £). This is reminiscent of the peak in mean-field calculations for 
the long-range spin-density wave ordered state. For fixed A (or doping), there is a spectral 
weight transfer from the high-frequency region to the low-frequency region as increasing e 
(or decreasing £). For the Hall conductivity lma xy (Q), we focused on the experimentally 
accessible frequency regime Q > 0.3t, and showed that Ima xy (Q) can be either positive or 
negative at small frequencies, depending on parameters A and Eq. A negative lma xy (Q) 
is rather non-trivial, and is a consequence of current vertex corrections.— 1 ^ For small Eq, 
liaa xy has a peak structure, reminiscent of the mean-field calculations. For fixed A (or 
doping), lma xy increases at low frequencies and decreases at high frequencies, as increasing 
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Eq (decreasing £, or increasing temperature). For fixed Eo, lma xy increases at low frequencies 
and decreases at high frequencies, as decreasing A (or increasing doping). 

In comparison to experiment, a xx calculated in our theory is ~ 2 — 3 times larger than 
data in the 0.1 < Q < 0.5eV range (Figure E]). We believe that this reflects the inadequate 
treatment of Mott correlations. Our calculated a xy is about 5 times larger than experiment 
(Figure [7] and Ref |48[). The structure, with a negative a xy at low frequencies and a positive 
value at higher frequencies is qualitatively consistent with data,— except that in our study, 
lm.a xy stays positive at high frequencies, unlike data.— Another minor point of difference is 
that in the data there is little temperature dependence of the zero crossing in a xy , while in 
the theory the zero-crossing point shifts with e - 

One advantage of the approach in the present paper is that the results for the electron 
spectral function, the conductivity, and the Hall conductivity are directly related to the 
mean-field results if the spin propagator takes the Kampf-Schrieffer form x ^ ^(q — Q)- 
However, this approach is insufficient for the study of transport properties in the low fre- 
quency limit. 
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Appendix A: Renormalization of the impurity relaxation time and the spin-fermion 
interaction vertex 

In our model, the fermions are scattered by both spin fluctuations, characterized by the 
interaction vertex g, and impurities, characterized by the relaxation time r. One important 
question is to study how one of the scattering process affects the other. This issue is 
addressed in this Appendix. We find that both renormalizations can be neglected in the 
sense to be discussed below. 
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(a) 



(b) 




FIG. 9: (a): the leading order diagram for the renormalization of the spin-fermion interaction 
vertex g due to impurity scattering, (b): the leading order correction to the impurity relaxation 
time due to spin-fermion interaction. The dashed line is the impurity line, the solid line is the bare 
electron propagator Go, and the wavy line is the spin- fluctuation propagator. 

1. Renormalization of the spin-fermion interaction vertex 

In this subsection, we discuss the renormalization of the spin-fermion interaction vertex 
g in the presence of impurity scattering. The leading order correction is given by Fig [9] (a), 



where u 2 is the impurity potential, and we only consider the static limit for the spin- 
fluctuation propagator. For the momentum transfer |q — Q| ~ £ -1 <C a~ l , the momentum 



pendently. As a result, Sg/g oc u 2 , and its dependence on ie n and q is estimated to be 
0(T/V, a|q — Q|) which can be neglected, where T> is the cut-off energy of the order of the 
fermion bandwidth. Applying the same argument to the impurity ladder diagrams, we find 
that the summation of the ladder diagrams gives a finite constant renormalization to g. As 
a result, we can neglect the diagrams that renormalizes g by properly redefining g. 

2. Renormalization of the impurity relaxation time r 

We now discuss the renormalization of the impurity scattering relaxation time r by the 
spin-fermion interaction. Figure |9] (b) shows the leading order term in calculating this 
renormalization, 




(Al) 



integral can be transformed to J de p de p+ Q, and the two integrals can be performed inde- 




(A2) 
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where 

T(q) = u 2 j T^CoCp', ze n ) 2 G (p' + q, ie n ). (A3) 

For the momentum transfer |q — Q| ~ <C a -1 , the J dp' integral can be transformed 
to J de p 'de p i + Q. Since the integral J de p > has a double pole, this leading renormalization is 
negligible. This argument persists to diagrams with more spin-fluctuation lines. Thus, the 
renormalization of the impurity scattering relaxation time due to the spin-fermion interaction 
is negligible. 



Appendix B: The static approximation to Eq ([6]) 

In this Appendix, we discuss the condition under which the static spin-fluctuation prop- 
agator Eq ([7]) can be used. For simplicity, we consider the electron self-energy in the leading 
order perturbation theory, 

E(p,je n )= 5 2 TV f -£L x (q iu n )- — . (Bl) 

' J {2ir) 2 ie n + iu n - e p+q 

ItOn I 

Substituting the spectral decomposition (see e.g. Ref [341]). 



dxlmx(q,x) 

X{q,iu n )= I : , (B2) 

j 7T X — lOJ n 



where x(q, u) is given by Eq ([6]), and summing over iu n , we obtain 



2 I dq f°° dx coth ^ - tanh e -^- 



E(p,ie B ) = ^ / j^- 2 / -^Im X (q,:r)— ^ (B3) 

y (2vr) 2 J„ 00 2?r ie n + x- e p+q 

The integral over x is restricted by Imx(q, x) to the region x < u; s f£ _2 . For w s f£~ 2 <C T, 
coth ^ « 2T/x, tanh can be neglected, and Eq (1B3I) is approximated as 

E(p, ie n ) « ^ 2 Txo / 777% f —JT—TTt — ~ r>^2^2 i — 2 — F~2 — A ' ( B4 ) 

J (2tt) 2 y 7T (1 + £ 2 (q - Q) 2 ) 2 + y 2 uj s £~ 2 y + ie n - 

where y = C, 2 x/u s f. Performing the y-integral by closing the contour to avoid the pole from 

the electron propagator, 

S(p, Z6„) = g 2 T J ( J )2 p + ( q _ Q)2 l£n + ^-2^ _ , p+g ' ( B5 ) 
where a = (1 + £~ 2 (q — Q) 2 ) is a quantity of order 1, and ia; s f£ _2 asgne n can be neglected 
compared to ie n . This leads to Eq (EJ) which was obtained in the static limit using Eq (J7J). 

For the current vertex function in the leading order perturbation theory, the use of Eq 
(J7J) is also justified in the same way as above; we can write an equation analogous to Eq 
f )28|) . and then split the two fermion Green's functions as was done in Sec IIII1 
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Appendix C: Summary of the formulas in the mean-field theory of the spin density 
wave state 



In this Appendix, we summarize the formulas of calculating the longitudinal and Hall 
conductivities in the mean-field SD W state. In calculating these quantities, the spin index a 
is irrelevant, giving an overall factor of 2, and will be neglected. The mean-field Hamiltonian 
is 

P V P 

with Q = (vr,7r), the two-component spinor = (c^c^q), and ii p = [ P ] . The 

\ A £p+qJ 

summation in the second equality is over the magnetic Brillouin zone as indicated by the 
prime. 

The imaginary-time (f ) electron Green's function in the mean-field theory is defined as 

G(p, f) ab = -< r T * p , a (f)¥j )6 (0) >, (C2) 
with the corresponding retarded function 

e - e p+Q + i/2r A 

A e - e p + i/2r / 

° (P ' e) = " _ u v + l/2T = ( e -E p + i/2T){e-E p+Q + i/2T)-^ (C3) 

where we have introduced a finite lifetime r. 

The real part of the longitudinal conductivity a xx (Q) is given by 

4*qY! £ ~ f Q iU + ^ Tr{^Im[^(p,c)]^Im[^(p,c + SI)]}, 

(C4) 

v x 

where = I p . Fig [10] (a) shows Kea xx (il) in the mean-field theory for A = 0.3t 



v p + Q, 



(solid line) and 0.6t (dashed line). 

The imaginary part of the Hall conductivity is given by^^i 

p s=l s=l 

(C5) 
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FIG. 10: The conductivities in the mean-field theory, (a): Hea xx (£l) obtained from Eq ()C4p with 
A = 0.3t (solid line) and A = 0.6t (dashed line), (b): Ima X y(Q) obtained from Eq (IC5h with 
A = 0.3t (solid line) and A = 0.6i (dashed line). 



where 



ImlT 



intra(l) 
R 



p,fl) 



imn 



intra(2) 
R 1 



p,fl) 



imn 



inter(l) 
R 1 



p,0) 



imn 



inter(2) 
R 1 



' — oo 



2^ 



[f(u)-f(u + Q)\ 

x |ReG + (p, u + n)A + (p, u + n)A + (p 

- ReG+(p, u) A + (p, u) A + ( Pl uj + fi) | , 

^ r ^im-f(u+Q)] 

J-oo 27T 

x {ReG~(p,tu + tt)A~(p,tu + tt)A~(p 

- ReG"(p, u)A~(p, u)A~(p, cu + ty 



-oo 27r 



[f(oj) - f(oj + n)] 



|ReG + (p, w + fi)A + (p, w + fi)A-(p, 

- ReG + (p,uj)A + (p,uj)A~(p,uj + 

/OO J 
_[/(«) -/(« + n)] 

[ReG"(p,a; + fi)A~(p,w + fi)A + (p, 

- ReG~(p, u)A~(p,uj)A + (p,uj + 



n)}, 



n)}, 



(C6) 



(C7) 



(C8) 



(C9) 
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imn 



inter(3) 
R 



2^ 



[f(u)-f(u + Q)]/2 



+ ReG~(p,u + fi)A + (p,w + fi) 
- (ReG + (p,u)A-(p,u) 



p,u) 



^mtra(l) = + ( E +x^ E +VV _ 2E+ x E+ y E+ xy 

^intra(2) = ( £ -»)2 £ ,- M + E'VV _ 2E~ X E~ y E~ Xy 



+ sin 2 26> p 



+ 



sin 3 20. 



a;\2 

pyp / > 



sin 2 2fl„ , kp . .. ((/^) 2 /if + (^) 2 /i™ - 2h;h y p h x p y ) 



sin 3 20 



A 



^inte r{ 3) = ^2 ^ ^ h Vf g ^ + (^f g yy _ 2 /£ 

+ sin 2 20 p + ft*^ - fc*^ - hfflh?) 



+ sin 2 26> p 



2/i. 



{{KfK x + (^) 2 ^ - *KK h ?) 



(CIO) 



+ ReG"(p, u)A + (p, w)J A + (p, w + n) 

/OO J 
^|/( w )-/(a; + n)]/2 

x |(^ReG + (p,a; + fi)A"(p,cj + fi) 

+ ReG"(p, w + n)A + (p, uj + 0)) A - (p, u) ( Cn ) 
- (ReG + (p,a;)A~(p,a;) 

+ ReG~ (p, cj) A + (p, w)) A~ (p, w + fi) I , 

(C12) 
(C13) 

(C14) 



(C15) 



(C16) 
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+ sin 2 26 p (hy p h; x + h x p g;h™ - h; g y p h;y - hy p h; y ) 

- sm 2 29 p J^ A2 ((W + (^) 2 < - iKK h 7) (C17) 

In the above, we have introduced notations G ± (p,e) = l/(e — £^ + i/2r) with = 
9 P ± v^p + A2 ' 9 P = (£ P + %+q)/2, h p = (e p - e p+Q )/2, tan# p = (h p - yfh% + A 2 )/A, 
^(pjO;) = — 2ImG ± (p, w), and have used the short-hand notation in which the superscript 
x, ?/ on the energy functions denotes derivative with respect to the corresponding momentum, 
e.g., E+ x = dE+/dp x , E+ xx = d 2 E+/dp 2 x) g x p = dg p /dp x , gf = d 2 g p /dp x dp yj Figure M 
(b) shows lma xy (Q) in the mean-field theory for A = 0.3t (solid line) and A = 0.6t (dashed 
line) . 

Appendix D: Explicit expression for Ima xy in Sec IIVI 

We now present the explicit expression for lma xy (Q) in Sec \1V\ 

lma xy (Q) = |*g^~ J ^f^{S? - ^ + Si + S* - S R }, (Dl) 



where 



with 



r°° (j f 

1mS?(p,u) = J —[f(e)-f(e + cu)]Re{U a -U i2 }, (D2) 



U n = [d y T x (p, e+ e + u> + )v* - d y T y (p, e+ e + cu + )v x p ] 

x [G R (p, e)d x G R (p, e + u)- G R (p, e + u)d x G R {p, e)] , (D3) 

t/12 = [^r x (p, e" e + u + )v* - d v T y {p, e~, e + u + )v x p ] 

x [G A (p, e)d x G R (p, e + u)- G R (p, e + u)d x G A (p, e)] , (D4) 

C/21 = [d x T x ( Pl e+ e + - ^T,(p, e+ e + u + )vf\ 

x [G R (p, e)0„G fl (p, e + to) - G R (p, e + co)d y G R (p, e)\ , (D5) 
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U 22 = [d x V x {p, e- e + - 9 x .r,(p, e", e + co + )v x p ] 

x [G A (p, e)a y G i? (p, e + u) - G R (p, e + co)d y G A (p, e)] , (D6) 

C/ 3 i = [r,(p, e + , e + - r tf (p, e + , e + 

x e + w)3,G*(p, e) - d x G R (p, e)d y G R (p, e + u>)], (D7) 

t/ 32 = [r x (p, e" e + u + )v y p - r„(p, e" e + w + )t£| 

x [9 x G fl (p, e + u;)^G A (p, e) - <9 x G A (p, e)d y G R (p, e + cu)], (D8) 

^41 = vy(e;i Q v y p+Q - e;i Q v; +Q )G R (p, e)G R (p, e + w )/(p, e + , e + (D9) 

^42 = <(^+Q< + Q - ^ Q < +Q )G' A (p, e)G*( P , 6 + W ) J( P , 6", 6 + (D10) 

tfo = - « +0 )G fi ( P , 6)G R (p, 6 + W ) J(p, e+ e + (Dll) 

^52 = <(^Q<+Q - 4VW GA (P' e ) G *(P> e + ^) J (P, e~, e + u + ). (D12) 
In writing these equations, we have used the notations e* 1 = e ± z5, and e + u; + = e + u + z5. 
The analytically continued vertex functions are 



T Q (p, e+, e + u + )= + ^(S«(p, e) - E«(p, e + w)) , (D13) 



and 



r Q (p, 6", e + w+) = < + (s A (p, e) - £*(p, e + w)) . (D14) 

The functions J(p, e^, e + w + ) are defined as 

f (P, «+, e + 0,+) = -2ffe €) + ^ ' + "> - 2 S " (P ' e) ~ E > ' + "> , (D15) 

and 

„ - , l A (p,e) +1 R (p,e + uj) o E A (p,e)-^(p,e + a;) 

J(p,e,e + uT) = — 2 — , D16 

uj + i/T (u + i/r) 2 

where 7 A,i? (p, e) = <9£ A,R (p, e)/de. The real part of cr-^fi) is obtained via the Kramers- 
Kronig relation, Eq (1T4"]) . with replaced by <t x?/ . 



1 A. W. Overhauser, Phys. Rev. 128, 1437 (1962). 



29 



2 N. P. Armitage, P. Fournier, and R. L. Greene, Rev. Mod. Phys. 82, 2421 (2010). 

3 A. Zimmers, J. M. Tomczak, R. P. S. M. Lobo, N. Bontemps, C. P. Hill, M. C. Barr, Y. Dagan, 
R. L. Greene, A. J. Millis, and C. C. Homes, EPL (Europhysics Letters) 70, 225 (2005). 

4 N. Bontemps, R. P. S. M. Lobo, A. F. Santander-syro, and A. Zimmers, Annals of Physics 321, 
1547 (2006). 

5 N. Bontemps, Physica C 460, 162 (2007). 

6 Y. Dagan, M. M. Qazilbash, C. P. Hill, V. N. Kulkarni, and R. L. Greene, Phys. Rev. Lett. 92, 
167001 (2004). 

7 S. Charpentier, G. Roberge, S. Godin-Proulx, X. Bechamp-Laganiere, K. D. Truong, 
P. Fournier, and P. Rauwel, Physical Review B 81, 104509 (2010). 

8 G. S. Jenkins, D. C. Schmadel, P. L. Bach, R. L. Greene, X. Bechamp-Laganiere, G. Roberge, 
P. Fournier, and H. D. Drew, Physical Review B 79, 224525 (2009). 

9 E. M. Motoyama, G. Yu, I. M. Vishik, O. P. Vajk, P. K. Mang, and M. Greven, Nature 445, 
186 (2007). 

10 Y. M. Vilk and A. M. S. Tremblay, Journal of Physics and Chemistry of Solids 56, 1769 (1995). 

11 Y. M. Vilk and A. M. S. Tremblay, Europhysics Letters 33, 159 (1996). 

12 P. A. Lee, T. M. Rice, and P. W. Anderson, Phys. Rev. Lett. 31, 462 (1973). 

13 M. V. Sadovskii, Soviet Physics - JETP 50, 989 (1979). 

14 J. Schmalian, D. Pines, and B. Stojkovic, Phys. Rev. B 60, 667 (1999). 

15 O. Tchernyshyov, Phys. Rev. B 59, 1358 (1999). 

16 M. V. Sadovskii and N. A. Strigina, Journal of Experimental and Theoretical Physics 95, 526 
(2002). 

17 G. Baym and L. P. Kadanoff, Phys. Rev. 124, 287 (1961). 

18 G. Baym, Phys. Rev. 127, 1391 (1962). 

19 O. K. Andersen, A. I. Liechtenstein, O. Jepsen, and F. Paulsen, Journal of Physics and Chem- 
istry of Solids 56, 1573 (1995). 

20 T. Moriya, Spin Fluctuations in Itinerant Electron Magnetism (Springer- Verlag, 1985). 

21 J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 

22 A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

23 L. B. Ioffe and A. J. Millis, Physics-Uspekhi 41, 595 (1998). 

24 A. Abanov, A. Chubukov, and J. Schmalian, Advances in Physics 52, 119 (2003). 

30 



25 I. Paul, C. Pepin, B. N. Narozhny, and D. L. Maslov, Phys. Rev. Lett. 95, 017206 (2005). 

26 H. v. Lohneysen, A. Rosch, M. Vojta, and P. Wolfle, Reviews of Modern Physics 79, 1015 
(pages 61) (2007). 

27 B. R. Patton and L. J. Sham, Phys. Rev. Lett. 31, 631 (1973). 

28 S. Takada and E. Sakai, Progress of Theoretical Physics 59, 1802 (1978). 

29 A. Larkin and A. Varlamov, Theory of fluctuations in superconductors (Clarendon Press, 2005). 

30 C. Caroli and K. Maki, Phys. Rev. 159, 306 (1967). 

31 G. Zala, B. N. Narozhny, and I. L. Aleiner, Phys. Rev. B 64, 214204 (2001). 

32 I. Paul, Phys. Rev. B 77, 224418 (2008). 

33 G. D. Mahan, Many-Particle Physics (Plenum Press, 1990), 2nd ed. 

34 A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinskii, Methods of Quantum Field Theory in 
Statistical Physics (Prentice Hall, 1963). 

35 B. L. Altshuler, L. B. Ioffe, and A. J. Millis, Phys. Rev. B 50, 14048 (1994). 

36 A. J. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 167 (1990). 

37 T. Moriya and K. Ueda, Advances in Physics 49, 555 (2000). 

38 S. Sachdev, Quantum Phase Transitions (Cambridge University Press, 1999). 

39 A. P. Kampf and J. R. Schrieffer, Phys. Rev. B 42, 7967 (1990). 

40 A. V. Chubukov and D. K. Morr, Physics Reports 288, 355 (1997). 

41 J. R. Schrieffer, Theory of superconductivity (Perseus Books, 1999). 

42 Z. Koba, N. Mugibayashi, and S. Nakai, Progress of Theoretical Physics 6, 322 (1951). 

43 H. Monien, Phys. Rev. Lett. 87, 126402 (2001). 

44 H. Fukuyama, H. Ebisawa, and Y. Wada, Progress of Theoretical Physics 42, 494 (1969). 

45 H. Fukuyama, Progress of Theoretical Physics 42, 1284 (1969). 

46 H. Kohno and K. Yamada, Progress of Theoretical Physics 80, 623 (1988). 

47 B. L. Altshuler and A. G. Aronov, Electron- Electron Interactions in Disordered Systmes (North- 
Holland, 1985), chap. Electron-electron interactions in disordered systems, pp. 1-154, Modern 
Problems in Condensed Matter Sciences, ISBN 0-444-86916-6. 

48 A. Zimmers, L. Shi, D. C. Schmadel, W. M. Fisher, R. L. Greene, H. D. Drew, M. Houseknecht, 
G. Acbas, M.-H. Kim, M.-H. Yang, et al., Physical Review B 76, 064515 (2007). 

49 H. Kontani, K. Kanki, and K. Ueda, Phys. Rev. B 59, 14723 (1999). 

50 G. S. Jenkins, D. C. Schmadel, P. L. Bach, R. L. Greene, X. Bechamp-Laganiere, G. Roberge, 

31 



P. Fournier, H. Kontani, and H. D. Drew, Phys. Rev. B 81, 024508 (2010). 
P. Voruganti, A. Golubentsev, and S. John, Phys. Rev. B 45, 13945 (1992). 



32 



